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ABSTRACT 

The destabilising effect of Hall diffusion in a weakly-ionised, Keplerian disc allows the 
magnetorotational instability (MRl) to occur for much lower ionisation levels than 
would otherwise be possible. However, simulations incorporating Hall and Ohm diffu- 
sion give the impression that the consequences of this for the non-linear saturated state 
are not as significant as suggested by the linear instability. Close inspection reveals that 
this is not actually the case as the simulations have not yet probed the Hall-dominated 
regime. Here we revisit the effect of Hall diffusion on the magnetorotational instability 
and the implications for the extent of MHD turbulence in protoplanetary discs, where 
Hall diffusion dominates over a large range of radii. 

We conduct a local, linear analysis of the instability for a vertical, weak mag- 
netic field subject to axisymmetric perturbations with a purely vertical wave vector. 
In contrast to previous analyses, we express the departure from ideal MHD in terms 
of Hall and Pedersen diffusivities jyn, and ryp, which provide transparent notation 
that is directly connected to the induction equation. This allows us to present a crisp 
overview of the dependence of the instability on magnetic diffusivity. We present an- 
alytic expressions and contours in the r]ii-r]p plane for the maximum growth rate and 
corresponding wave number, the upper cut-off for unstable wave numbers, and the 
loci that divide the plane into regions of different characteristic behaviour. We find 
that for sign(i?z)?7H < — 2w^/r2, where va is the Alfven speeds and ft is the Keplerian 
frequency. Hall diffusion suppresses the MRI irrespective of the value of rfp. 

In the highly-diffusive limit the magnetic field decouples from the fluid perturba- 
tions and simply diffuses in the background Keplerian shear flow. The diffusive MRI 
reduces to a diffusive plane-parallel shear instability with effective shear rate We 
give simple analytic expressions for the growth rate and wave number of the most 
unstable mode. 

We review the varied and confusing parameterisations of magnetic diffusion in 
discs that have appeared in the literature, and confirm that simulations examining 
the saturation of the instability under Hall-Ohm diffusion are consistent with the 
linear analysis and have yet to probe the "deep" Hall regime |?7h| > r]p > v\/n 
characteristic of protoplanetary discs where Hall diffusion is expected to overcome 
resistive damping. 

Finally, we illustrate the critical effect of Hall diffusion on the extent of dead 
zones in protoplanetary discs by applying a local stability criterion to a simple model 
of the minimum-mass solar nebula at 1 au, including x-ray and cosmic-ray ionisation 
and a population of 1 /im grains. Hall diffusion increases or decreases the MRI-active 
column density by an order of magnitude or more, depending on whether B is parallel 
or antiparallel to the rotation axis, respectively. We conclude that existing estimates 
of the depth of magnetically active layers in protoplanetary discs based on damping 
by Ohm diffusion are likely to be wildly inaccurate. 

Key words: accretion, accretion discs - instabilities - MHD - planetary systems: 
protoplanetary discs - stars: formation 
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1 INTRODUCTION 

The mechanism by which angular momentum is lost by the 
material in accretion discs is a classic problem in astro- 
physics. While gravitational torques associated with a bi- 



nary companion (e.g. Sawada et al. 19871 or gravitational 
instabilities in massive discs ( [Adams et al.||1989| may play 
an occasional role, magnetic tension appears likely to be the 
most common source of the torque required to remove an- 
gular momentum from the material in the disc and permit 
it to spiral towards the central object. 

Magnetically-mediated accretion disc models fall into 
two broad paradigms. In disc-jet models, a large-scale 
poloidal field threads the disc and the magnetic tension of 
the field lines anchored to the disc accelerates material out- 
ward and upward along the field lines away from the disc 



surfaces to form a collimated bipolar jet ( Blandford & Payne 
|1982[ ). The net effect is to transfer angular momentum from 
the disc material to the accelerated outflow. In the absence 
of an outflow the disc may still be braked through torsional 
magnetic coupling to an extended envelope or surrounding 
cloud (e.g. Krasnopolsky fc Konigl|2002 L 

In models invoking magnetically-driven turbulence, on 
the other hand, tension in a weak magnetic field viru- 
lently destabilises the orbital shear flow in the disc via the 



magnetorotational instability (MRI; Velikhov 1959 Chan- 
drasekhar|1960| [Balbus fc Hawley|1991| ), driving MHD tur 
bulence that transports angular momentum radially out- 
wards (e.g. |Balbus fc Hawlcy 1992} [Stone et al.|1996[ ). Gen- 
erally, one expects that the jet models apply when the mag- 
netic and thermal pressures are comparable at the disc mid- 
plane, while substantially sub-thermal magnetic fields are 
subject to the MRI and are incapable of driving strong jets. 
It is conceivable that these two transport mechanisms may 
operate at the same radius for intermediate strength fields, 
with the MRI operating for a range of z adjacent to the 
mid-plane of the disc and vertical transport dominating at 
greater heights (e.g. |Salmeron et al.|2007| . Similarly, as the 
field strength is expected to decrease with radius, there could 
also be an inner region where jets dominate, with the MRI 
operating at larger radii ( [Combet fc Ferreira|2008| ). 

In this paper we focus on the efficacy of the MRI- 
driven turbulence in protostellar/protoplanetary discs. In 
this context, apart from transporting angular momentum, 
MRI-driven turbulence selectively heats the disc, effectively 
mixes chemical species (e.g. Semenov et al.'2006'), and affects 
the transport, aggregation, and sedimentation of dust grains 



i Johansen fc Klahr|2005l [Turner et al.|2006|[Fromang fc Pa- 



paloizou[|2006 1, the assembly blocks of planetesimals in the 
'core accretion' scenario of planet formation. Furthermore, 
the effective viscosity of the disc associated with MHD tur- 
bulence determines the ability of protoplanets to open gaps 
in the disc and also contributes to the delicate imbalance 
of gravitational torques that determines the rate and di- 



rection of planetary migration e.g. (Matsumura fc Pudritz 
2003} [Nelson fc Papaloizou|2004[ [Johnson et al.|2006| |. 

It is of great importance, therefore, to determine the 
location, extent and behaviour of the magnetically-active 
regions of protoplanetary discs, which are restricted by the 
very low fractional ionisation beyond ~ 0.1 au from the cen- 
tral star, where thermal ionisation is ineffective ( |Gammie| 
19961. The dominant sources of ionisation - X-rays or UV 



radiation from the forming star, and possibly interstellar 
cosmic rays - are unable to penetrate to the disc mid-plane 
except perhaps at larger radii where the column density is 
low. At intermediate radii a very low level of ionisation is 
maintained by cosmic rays or, if these are absent, the decay 



of radioactive elements mixed with the gas (Semenov et al. 
2004| [Glassgold et al.[[2005| ). MRI-driven MHD turbulence 



is therefore thought to be restricted to the surface layers 
(IG ammie Il996| ) over the range ~ 0.3-20 au where shielding 



of external ionisation sources results in very low levels of 
ionisation. 

Assessments of the location and extent of magnetically- 
active zones (or equivalently, their complementary magnetic 
'dead zones') typically invoke only damping of the linear 



MRI by Ohm diffusion (e.g. 


Hayashi 1981! |Jin 


1996 


Fro- 


mang et al.|2002| [Matsumura 


fc Pudritz.2003|[Ilg 


ner fc Nel- 


son[[2006| Turner fc Sano[[2008|. However, Ohm diffusion is 



the dominant magnetic diffusivity only at the very high den- 
sities achieved close to the disc mid-plane within 1 au of the 
central star. Instead, Hall diffusioij^ dominates at the inter 



mediate mid-plane densities between 1 and 30 au (Wardle 
fc Ng[[1999||Sano fc Stone|[2002a[ ), and ambipolar diffusion 



beyond about 30 au. Note, however, that as the density de- 
creases towards the disc surfaces, different diffusivities may 
dominate in different layers. The neglect of ambipolar diffu- 
sion in assessing magnetic activity is less drastic, as it domi- 
nates near the disc surface where the fractional ionisation is 
relatively high and magnetic diffusivity may be insufficient 
to stabilise the MRI (e.g. |Perez-Becker fc Chiang[[2011 K 

Hall diffusion provides a dissipation-free pathway for 
the magnetic field that either enables the MRI to grow 
despite the damping effect of Ohm or ambipolar diffusion 
or suppresses it entirely (Wardle 19991 Balbus fc Terquem 



20011, depending on the orientation of the magnetic field. 
Based on the linear growth rates Hall diffusion may dras- 
tically extend or restrict the reach of magnetic activity in 



protoplanetary discs (Wardle 20071 and hkely modifies the 



transport and dissipative properties of the resulting turbu- 
lence. However, pioneering shearing-box calculations includ- 
ing both Ohm and Hall diffusion did not appear to confirm 
this expectation, suggesting instead that Hall diffusion was 



unable to counter the damping effect of Ohm diffusion ( Sano 
fc Stone|2002a|b I. As a result, subsequent evaluations of the 



extent of dead zones in protoplanetary discs are based either 



on the damping effects of Ohm diffusion (e.g. Turner fc Sano 
2008[ [Turner fc Drake|2009[ [Kretke fc Lin|2010[), or on O hm 
and ambipolar diffusion ( [Perez-Becker fc Chiang[[20lT| [Bai| 
20TT|). 



The apparent failure of fully-developed turbulence to re- 
spect the underlying linear instability in the Hall-diffusion- 
dominated limit is surprising, as it does so in the Ohm 



(e.g. Turner fc Sano||2008j and ambipolar diffusion regimes 
Stone 2011|. However, closer inspection of the 



(Bai 



parametrisation of the diffusive effects adopted by [Sano fc[ 
Stone (2002a b) reveals that their simulations did not probe 



the important regime in which (a) Hall diffusion dominates 
Ohm diffusion and (b) both dominate inductive effects. In- 



^ Despite the lack of dissipation associated with Hall drift we 
adopt the term "Hall diffusion" to maintain consistency with the 
notion of a tensor diffusivity appearing in the induction equation. 
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deed, none of the results of Sano & Stone ( 2002a|b' I conflict 
with expectations based on the linear analysis. The confu- 
sion that has arisen can largely be attributed to the adoption 
of a magnetic Reynolds number Rcm and a Hall parameter 
X to characterise the magnitudes of Ohm and Hall diffusion 
respectively. Then the criteria for Ohm and Hall diffusion 
to dominate the inductive term are Rcm < 1 and X > 2, 
respectively. Crucially, the ratio of Hall to Ohm diffusion is 
then ^XRcm- This is never ^ 1 when Rcm < 1 in the sim- 
ulations, apart from the zero-net-flux simulations S14-S16 



of Sano & Stone ( 2002b I in which Hall diffusion suppresses 



the MRI in the Bz < regions, as one might expect based 
on the linear analysis. 

The purpose of this paper is to emphasise the impor- 
tance of Hall diffusion for MRI-driven turbulence in proto- 
planetary discs. We do this by considering the simplest ex- 
ample of the MRI - perturbations with wave vectors parallel 
to an initially vertical magnetic field - and demonstrating 
that the column density of the MRI-unstable region at 1 au 
in the minimum mass solar nebula increases or decreases 
(depending on the sign of Bz) by an order of magnitude 
when Hall diffusion is included (see Fig. |14[ ). We begin by 
describing the field-line drifts induced by ambipolar, Hall, 
and Ohm diffusion in §2 and give a qualitative discussion 
of their effect on the MRI. In §3 we present an overview of 
the results of a linear analysis of the MRI in the presence 
of diffusion and give a coherent survey of its dependence on 
Pedersen (ie Ohm-|-ambipolar) and Hall diffusion by consid- 
ering contours of the growth rate and wave number of the 
fastest growing mode and the range of MRI-unstable wave 
numbers in §4. We then compare previous notation with the 
diffusivity notation we now favour and discuss the numeri- 
cal simulations of Sano & Stone ( 2002a|b |. In §5 we estimate 
the column density of the MRI-unstable layer at 1 au in the 
MMSN using the criterion that there should exist a local un- 
stable MRI mode with wave number satisfying kh > 1. The 
results are dramatic, indicating that Hall diffusion increases 
the column density of the MRI-unstable surface layers by 
an order of magnitude depending on whether Bz is positive 
or negative, respectively. Finally, we summarise our results 
and conclusions in 56. 



2 MAGNETIC DIFFUSION AND THE 

MAGNETOROTATIONAL INSTABILITY 

2.1 Magnetic diffusion 

The magnetic field evolves according to the induction equa- 
tion 



dB 

dt 



V X (v X B) - cV X E' . 



(1) 



where E' is the electric field in the local instantaneous rest 
frame of the fiuid, which satisfies the generalised Ohm's Law 

cE' = 77a(V X B)x + ??h(V X B) X B + rj{V X B) , (2) 

where we use subscripts || and _L to refer to the orientation 
with respect to the local magnetic field, and 77A, Vn, and 77 
are the ambipolar. Hall, and Ohm diffusivities, respectively. 
It will prove useful to recast the induction equation in a form 
that makes explicit the drift of the magnetic field through 



the fluid, i.e 

dB 

~dt 
where 

E' X B 



V X [{v + vb)x B - ?7(V X B) 



(3) 



52 



(V X B)x X B (V X B) 

= ^ m — ^ 

and the Pedersen diffusivity 77P is given by 
?7p = ?7A + »7, 



(4) 



(5) 



This form makes explicit the role of the diffusivities in trans- 
porting flux, and emphasizes that only Ohm diffusion is ca- 
pable of destroying magnetic flux (i.e. via the 7;(V X B)| 
term) . 

The diffusivities are determined by the abundances of 
charged particles and their collision cross-sections with the 
neutrals (e.g. Cowling 1957 Wardle & Ng 1999 1, and - when 



the fractional ionisation is not too small - with each other. In 
general, r), rjA, and t^p are all positive, whereas rju is positive 
or negative depending on whether positively- or negatively- 
charged species, respectively, are on average more decoupled 
from the magnetic field by collisions. For example, for a 
weakly-ionised ion-electron-neutral plasma we obtain 



c rae"iep 
iTre^ne 



and 



B' 
A-K-yi ppi 



cB 



(6) 



(7) 



(8) 



Air e He 

(e.g |Konigl|1989[|Balbus fc Terquem|2001| ). The diffusivities 
exhibit more complex dependencies on density, field strength 
and electron abundance when charged grains are present. In 
this case the Hall diffusivity is negative at relatively low 
densities when neutral collisions decouple grains from the 
magnetic field, but not the electrons and ions, and is positive 
at higher densities when ions are also decoupled from the 



magnetic field (e.g. Wardle & Ng 19991. The diffusivities 
exhibit significant spatial gradients if the ionisation rate and 
gas density do so, as in protoplanetary discs (e.g. [Wardle 



[2007'). Fortunately, for a linear analysis the response of the 
diffusivities to perturbations is not needed as V X B vanishes 
in the initial equilibrium state. Thus in our local analysis 
we can regard the diffusivities as specified constants, but 
we shall return to them when considering the application to 
protoplanetary discs in !|5] 

2.2 Effect on the MRI 

Before deriving the dispersion relation for the magnetorota- 
tional instability in diffusive media we first give an intuitive, 
physical description of the effect of magnetic diffusion on its 
growth and properties. 

We first examine the growth of the instability in a Ke- 
plerian accretion disc under ideal-MHD conditions and as- 
suming an initially vertical magnetic field (see Fig. [l] top- 
left panel). We consider, in particular, the situation where 
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JxB 




Bq+6B 



Figure 1. A sketch of the development of the magnetorotational 
instabiUty (MRI) in an initially vertical magnetic field Bq embed- 
ded in an unstratified, Keplerian disc for a vertical wave vector. 
Top panel shows a poloidal view of (left) the initial field con- 
figuration and (right) the perturbed field and associated current 
density (J), Lorentz force (J X B) and deviation from local Kep- 
lerian velocity. The lower panel shows in plan view the perturbed 
fluid elements at points 1 and 2. At point 1 the magnetic ten- 
sion provides some radial support against gravity and supplies a 
negative torque. As a result the material here is sub- Keplerian 
and moving inwards, as indicated by the green vector which indi- 
cates the departure from Keplerian rotation. The field transfers 
the angular momentum to the fluid at locations such as 2 where 
the field is buckled outwards, and the fiuid there spirals out. The 
differential radial motion of the fluid that is losing and gaining 
angular momentum at 1 and 2 respectively enhances the buckling, 
leading to runaway growth. 



alternate layers of fluid are displaced from their equilibrium 
position, such that point 1 moves radially inwards (and for- 
wards in azimuth) while point 2 is shifted radially outwards 
(and backwards in azimuth) . In this case the magnetic field is 
buckled because the field lines are frozen into the fluid (top- 
right panel). The deformation of the lines creates magnetic 
tension forces which are directed outwards and backwards in 
azimuth (at point 1), or inwards and forwards in azimuth (at 
point 2), as depicted in the lower panel of Fig. 1. At point 
1 then the magnetic tension provides some radial support 
against gravity and supplies a negative torque. This means 
that the fluid element at point 1 will lose angular momen- 
tum and spiral inwards. The opposite takes place at point 2, 
with the fluid element there spiralling outwards. As a result, 
the magnetic field buckling and associated tension increases, 
leading to runaway growth. The fastest growing mode has 
growth rate | f2 on scales with wave numbers k ~ fl/vA, 
where Q — VK/r is the Keplerian frequency and va is the 
Alfven speed ( |Balbus fc Hawley|199l"| . 

Magnetic diffusion modifies the above picture by allow- 
ing slippage between the field lines and the fiuid, so that 
there is no longer a direct connection between the relative 
displacement of the fluid layers and the buckling of the mag- 
netic fleld lines. One might expect that this effect should al- 
ways reduce the instability because the displacement of the 
field would lag that of the fiuid. However, this intuition is 
based on the limit of Ohm or ambipolar diffusion, in which 
the drift is in the direction of the magnetic stresses and tends 
to straighten up field lines. Hall diffusion, by contrast, cre- 
ates a drift orthogonal to the tension forces and may, there- 




V + V 



z*- r 




Figure 2. Plan view of the effect of ambipolar or Ohm diffusion 
on the development of the MRI. The current density and associ- 
ated magnetic stress on the fluid at points 1 and 2 of Fig. ^ are 
indicated in blue. Because of flnite diffusivity the fleld lines drift 
through the fluid at velocity hb {red vectors) in the direction 
parallel to J X B (see eq [ojl . The net drift of the fleld line with 
respect to Keplerian rotation, v + is indicated by the black 
vector. In this case the net drift is marginally inwards at point 
1. At point 2, the current density and magnetic stresses are re- 
versed, so the field drifts in the opposite direction. As a result the 
field lines drift inwards and outwards at reduced rates at points 1 
and 2 respectively, and the instability proceeds more slowly than 
in ideal MHD. 



fore, enhance or suppress the radial buckling depending on 
the situation. It is this feature that gives Hall diffusion its 
unique properties and that we aim to explain here. 

In the ambipolar diffusion limit, the magnetic field 
is frozen into the ions and electrons, which drift together 
through the neutral component of the fiuid. Collisions with 
the neutrals then transmit magnetic stresses to the bulk of 
the gas. As the fractional ionisation is low, the ion and elec- 
tron inertia and thermal pressure are negligible and the drift 
velocity of the ions, electrons and field is determined by the 
balance between the Lorentz force on the ions and the col- 
lisional drag with the neutrals, 

(9) 



Vp — Vi — V = 



and so is parallel to J X i?. In eq. (|9|, p and pi 
neutral and ion density, respectively, and 



rrii + m 



are the 



(10) 



where {av)i measures the rate coefficient of momentum ex- 
change via collisions with the neutrals, taken to have mean 
mass m, and nii is the ion mass. When J ■ B — 0, such 
as envisaged here. Ohm diffusion also produces a drift in 
the same direction, so we have adopted the subscript P to 
denote the Pedersen drift (see eq |A10 1. 

Fig.[2]shows the field line drift (red vectors) for the fluid 
conflguration depicted in Fig. [T] but now incorporating the 
effect of ambipolar or Ohm diffusion. The black vectors in- 
dicate the total drift of the magnetic field in the local Keple- 
rian frame, and correspond to the vector sum v + vb (see eq 
AlO I. In this case the net effect of the drift is to reduce the 



radial and azimuthal stretching of the field, and therefore - 
as one might guess - reduce the degree of instability. 
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V + v„ 




Figure 3. Plan view of the effect of magnetic diffusion [incorpo- 
rating the ambipolar, Ohm and Hall (> 0) terms] on the develop- 
ment of the MRI. The current density and associated magnetic 
stress on the fluid at points 1 and 2 of Fig. [l] are indicated in 
blue. At point 1 the magnetic tension provides some radial sup- 
port against gravity and supplies a negative torque. As a result, 
the material here is sub-Keplerian and moving inwards, as indi- 
cated by the green vector (which denotes departure from Kep- 
lerian rotation). Because of the finite diffusivity, the field lines 
drift through the fluid [red vectors) with a velocity determined 
by the magnitudes of the ambipolar. Ohm and Hall diffusivities 
(see eq |A10| . Ambipolar and/or Ohm diffusion contributes a drift 
vp parallel to J X B, while the Hall drift i^jj is antiparallel to 
the current density J. The net drift of the field line with respect 
to the Keplerian rotation, v -f v^, is indicated by the black vec- 
tor. In this case the net drift is inwards. At point 2, the current 
density and magnetic stresses are reversed, so the field drifts in 
the opposite direction. As a result the field lines drift inwards at 
point 1 and outwards at point 2 and the instability proceeds. 



In the Hall limit, collisions with the neutral gas are 
frequent and strong enough to decouple the ions from the 
field, but the electrons - which have a smaller collision cross 
section with the neutrals and a higher charge-to-mass ratio 
- remain well coupled. The magnetic field is then frozen into 
the electrons, and these drift together through the neutrals 
and ions, which remain tightly-coupled by collisions. In this 
limit, the field drift speed through the neutrals is given by 
the ion-electron drift and hence is antiparallel to the current 
density, 

J 



(11) 



In eq. rie (= n^) is the electron number density and e 

is the elementary electric charge. 

More generally, when all diffusion mechanisms are taken 
into account the drift velocity vb of the field lines through 
the fluid is the sum of these two orthogonal contributions 



(see eq AlO below). The implications for the magnetorota- 
tional instability are sketched in Fig. |3] which now includes 
a Hall drift antiparallel to J. The generic geometry of the 
MRI means that the Hall drift at point 1 is directed inwards 
and retrograde in azimuth, and it is outwards and prograde 
at point 2. The net effect is to exacerbate the radial buck- 
ling of the field, but to reduce the buckling in the azimuthal 
direction. We show in the next section that it is the radial 
effect that is critical, as azimuthal field is always created 





Figure 4. As for Fig. [3] but after a global reversal of the mag- 
netic field, as would be expected if the initial magnetic field was 
antiparallel to the disc's angular velocity vector. In this case, the 
magnetic stresses are unaffected — and so, therefore, are the fluid 
velocity and Ohm/ambipolar drift current density. However, the 
current density is reversed, and therefore so is the Hall drift. The 
magnetic field now drifts outwards at 1 and inwards at 2, reduc- 
ing the radial buckling. The MRI is partly or entirely suppressed 
in these circumstances (see text). 



out of the underlying Keplerian shear of the disc. So Hall 
diffusion in this case (e.g. when the magnetic field is aligned 
with the angular velocity vector of the disc) tends to be 
destabilising. 

This is not always the case, however, because the direc- 
tion of the current density, and therefore of the Hall diffu- 
sion, reverses upon global reversal of the magnetic field. This 
situation is illustrated in Fig. |4] The magnetic stresses, the 
fiuid velocity, and the field-line drift associated with Ohm 
and ambipolar diffusion are not affected by this reversal. 
Note, however, that Hall diffusion now acts to stabilise the 
disc by acting against the radial buckling of the field that 
would otherwise be driving the fluid motions. 

To summarise. Ohm and ambipolar diffusion are stabil- 
ising effects, whereas Hall diffusion may be destabilising or 
stabilising depending on whether the initial vertical mag- 
netic field is parallel or antiparallel to the rotation axis, re- 
spectively. This asymmetry refiects the fundamental asym- 
metry in the microscopic properties of the positive and neg- 
ative charged species in the fluid in this limit. 



3 LINEAR ANALYSIS OF THE MRI 

We consider a small region of an axisymmetric, geometri- 
cally thin and nearly Keplerian disc, threaded by a vertical 
magnetic fleld, with sound and Alfven speeds (cs and va) 
at the mid-plane that are both small compared to the local 
Keplerian speed vk- We assume that radial gradients are on 
the scale of r and neglect vertical stratiflcation of the initial 
equilibrium state, so that our analysis only holds near the 
mid-plane at heights z <C Ca/il. The initial state is in Keple- 
rian rotation with a uniform density, pressure, and vertical 
magnetic field B — sBz (where s = ±1). We linearise the 
equations around this state and seek solutions for axisym- 
metric perturbations of the form exp(i^f — ikz). The equa- 
tions for perturbations in density, pressure and form a 
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separate system that describes vertically-propagating sound 
waves. The system of linear equations in the remaining per- 
turbations involve fluctuations in the r and (j) components 
of B, V and vb- 

A derivation of the dispersion relation is outlined in Ap- 
pe ndices \K\ and [B| This dispersion relation was first derived 
by Wardle (19991, and later extended to more general ge- 
ometries in the Hall-Ohm limit ( Balbus &: Terquem||2001 |, 
and then including ambipolar diffusion (Desch 20041; the 



ambipolar diffusion limit was also considered by |Kunz fc| 
Balbus (20041. Our emphasis here is on the dependence of 



the instability on the Pedersen and Hall diffusivities 77P and 
?7h; a detailed analysis is presented in the appendices. An 
overview is provided by Fig. [5] which illustrates the qualita- 
tive changes in the growth rate vs wave number curves for 
different choices of the diffusivities. 

First, note that the ideal-MHD limit holds at the origin 
(i. e. 77H = ??p = 0), the Hall MHD limit holds along the 
horizontal axis (r;p = 0) and the Ohm or ambipolar diffusion 
limits hold along the vertical line (?7h = 0). Recall also that 
only the half plane r/p > is physically relevant. It turns 
out that there are three distinct forms of the resulting i'(fc) 
curve, corresponding to regions labelled I, H and HI in the 
plane, as illustrated in Fig. [S] In region I wave numbers less 
than a cutoff kc are unstable with the maximum growth rate 



attained at an intermediate wave number fco (see eqs B12 



B13 and B14I. In region II all wave numbers are unstable 



with the maximum growth rate still occurring at finite wave 
number. Finally, in region III all wave numbers are unstable 
and the maximum growth rate is approached asymptotically 
as fe — >■ 00. 

Having delineated these three regions, we now consider 
how the critical wave number kc, fastest growth rate 
and corresponding wave number ko vary across the entire 
?7p-'7h plane. Contours of the growth rate and wave num- 
ber of the most unstable mode are plotted in this plane in 
Fig. [5] The growth rate increases clockwise, from f2 along 
the vertical line srju ~ —2 up to 0.75 f2 for the horizontal 
line ?7p = for srju > —4/5. In the absence of Hall diffu- 
sion, the maximum growth rate vo declines with increasing 
(Ohm and/or ambipolar) diffusivity (e.g. moving vertically 
upwards through the s-qn = point in the horizontal axis), 
with vo ~ l^p^ for ^7? ^ 1. The most important effect of 
Hall diffusion, apparent from Fig|6] is that the growth rate 
of the MRI exceeds 0.3 fl for srju ^??p, even for arbitrarily 
large rjp. More generally, the addition of Hall diffusion at 
fixed rjp increases the growth rate if srju > and decreases 
it when sr/n < 0. For large values of rjp, eq ( |B13[ | shows that 
srju /rjp « 24i/o/(9 — 16uq). It is this fact that has the po- 
tential to modify the extent of dead zones in protoplanetary 
discs, as we explore in ^ 

The wave number of the fastest growing mode (blue 
contours in Fig. [6| decreases as the diffusivity is increased. 
Again, the contours are not arranged so that the highest 
wave numbers occur in the ideal-MHD limit, but to the 
srjH < side, within the boundary between regions II and 
HI (traced by the fco = 00 contour). 

Contours of constant kc are semicircles, as plotted in 
Fig. [7] While the range of unstable wave numbers is reduced 
for large values of srju and rjp, as one might expect, the range 
is not maximised in the ideal limit (ie. at the origin) but in 
regions II and HI, bounded by the kc = 00 contour. 



In general we note that increasing rjp decreases the 
maximum growth rate and the characteristic wave numbers, 
whereas increasing srju above —2 increases the maximum 
growth rate and may either increase (when S77H + 1 ^ f?p) or 
decrease (when srju + ^^Vp) the corresponding wave num- 
ber. 

Overall, these patterns place the ideal. Ohm (or am- 
bipolar) and Hall regimes in context, and for the first time 
we see an overview of the effect of magnetic diffusivity on 
the linear MRI. In particular, there is nothing special about 
the Ohm/ambipolar limit, e.g. the behaviour of the instabil- 
ity in the presence of diffusion is not qualitatively different 
for srju = 2 vs ?7h =0. Even the ideal-MHD limit does not 
stand apart as remarkable, although it still holds a special 
place conceptually because fiux freezing holds and it is easier 
to visualise. 

What does stand out is the part of the plane in the lower 
left, regions II and III, characterised by high wave numbers 
and the spraying out of the growth contours. In this part of 
the plane, the instability operates in the "cyclotron limit" 
rj± ^ 1 and fc^ 3> 1 (see appendixjcjl , in which the instability 
arises in the competition between magnetic diffusion and 
advection of the field by the fiuid; generation of from 
Br by the Keplerian shear ffow is negligible. The lack of 
any k dependence in this regime occurs because both the 
magnetic diffusion and the magnetic stresses on the fluid 
(which are responsible for the fluid displacement) scale as 
fc^. This short- wavelength, low-frequency limit corresponds 
to the cyclotron mode of the magnetised fluid , which has 



frequency ujh = v^/lrjul (Wardle & Ng 



1999 



Pandey & 



Wardle 2008 1. This mode couples effectively to the Keplerian 



rotation as long as the sense of circular polarisation matches 
that of epicyclic motion, i.e. as long as B^rju < 0- The other 
short-wavelength mode, the high-frequency whistler {uj ~ 
k'^v\/ujH), is unable to couple effectively to the rotation 
( Wardle|1999t . 

The ma^ximum growth rate contours emerging from re- 
gions II and HI continue on to attain another limit when 
')p + ^ in which the field evolves in response to shear 
and diffusion, without significant feedback from the pertur- 
bations that it induces in the fluid flow. Instability in this 
case relies on the keplerian shear flow generating B^ from 
Br, and the tendency of Hall diffusion to convert B^ back 
into Br - This brings the potential destabilising effect of Hall 
diffusion in shear flows ( Kunz|[2008 1 to the fore and shows 
that it is quite independent of rotational effects - i.e. the 
Coriolis and centripetal acceleration - that drive the MRI. 
In Appendix [C] we obtain simple analytic expressions for 
growth rate in plane-parallel shear flows as a function of k 
for arbitrary diffusivity, the results are plotted in Fig. [8] For 
more general field and wave vector configurations, ambipolar 
diffusion plays a similar role, albeit hindered by dissipation 
(Kunz & Balbus||2004l |Desch||2004l |Kunz||2008l Pandey & 



Wardle, in prep). 



4 COMPARISON WITH NUMERICAL 
SIMULATIONS 

In this section we explore the relationship of our linear cal- 
culations to the nonlinear, unstratified shearing box simula- 
tions of MRI-driven turbulence by Sano & Stone ( 2002a|b I , 
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Figure 5. Schematic dependence of the behaviour of the MRI on the Hall and Pedersen (Ohm+ambipolar) diffusivities, rjn and ryp, for 
an initial magnetic field sBz in a Keplerian disc, where s = ±1 is the sign of B^. The field is assumed to be weak so that stratification 
can be neglected. Physical values of the Pedersen conductivity rjp are non-negative, and there are no unstable modes for Hall diffusivities 
sriiifl/v^ < —2. The unstable region srjufl/v'^ > —2 is subdivided according to the dependence of the growth rate on wave number k, 
which is sketched in the insets labelled I— IH. In region I (outside the red locus given by eg. |Bll[ l, wave numbers less than a cutoff fcc 
are unstable with the maximum growth rate i/q attained at fco (see eqs |B12| |B13[ and |B14] |. Between the red and blue loci (region II) 
all wave numbers are unstable with the maximum growth rate still occurring at finite wave number. Within the blue locus given by eq. 
|B15| (region III) all wave numbers are unstable and the maximum growth rate is approached asymptotically as fc — > oo. 




I I 1 I 1 I 1 I 

-2 2 4 

s riH ^ / 

Figure 6. Contours of the maximum growth rate of the MRI (solid black lines, units fl) and corresponding wave number (dashed 
blue lines, units Sl/ii^) as a function of Hall and Pedersen (Ohm+ambipolar) diffusivities. The innermost blue contour shows where 
fco becomes infinite, corresponding to the boundary between regions II and III in Fig. [5] Within this region the fastest growth rate is 
approached asymptotically as fc — > oo; this transition is responsible for the curvature of the black contours within this region (see eqs 
|B13| and [BT6l l. 
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s riH 5^ / V/ 



Figure 7. Contours of the critical wave wave number kc below 
wiiich instability sets in (units VI/va)- The red innermost contour 
shows where kc becomes infinite and all wave numbers are unsta- 
ble, corresponding to the boundary between regions I and II in 
Fig. [5] 
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Figure 8. Growth rate versus wave number for the Hall-driven 
diffusive instability for a mutually orthogonal magnetic field, 
shear velocity, and velocity gradient v' . For Keplerian rotation 
and cylindrical geometry, the effective velocity gradient is v' = 
The curves correspond to different values of the ratio of the 
Pedersen and Hall diffusivities. 



which explicitly included both Ohm and Hall diffusion. 
These authors' results could be interpreted as implying that 
the presence of Hall currents has little effect on the criti- 
cal (minimum) degree of magnetic coupling (the coupling 
between the neutral gas and the magnetic field) for the in- 
stability to operat^ We show, however, that their calcula- 
tions do not yet probe deeply enough into the Hall regime 
for Hall diffusion to be able to impact on the development 
and properties of the instability. In fact, their solutions are 
in agreement with expectations from local and stratified lin- 
ear analyses for the same values of the parameters describing 



^ This field-matter coupling parameter is often referred to as the 
magnetic Reynolds number Rem- As we expound below, however, 
it really corresponds to the Elsasser number A (see discussion 
below and Table 



the initial fluid conditions. Finally, we use the ratios of the 
different terms in the induction equation in order to delin- 
eate the region of parameter space where Hall diffusion is 
expected to substantially modify the growth rate, and spa- 
tial scale, of the MRI-unstable modes. First, however, and 
for the sake of clarity, we discuss the notation used in the 
literature to characterise the fluid, and its magnetic activity, 
in diffusive environments. 



4.1 Fluid parameters 

The growth and structure of the MRI are strongly depen- 
dent on the magnetic fleld strength, geometry, and the na- 
ture - and magnitude - of the magnetic diffusivity. Different 
parametrizations have been used in the literature to char- 
acterise these properties, depending on the formulation of 
the problem and the adopted form of the induction equa- 
tion. Some authors use a tensor diffusivity, and the induction 
equation then takes the form 

dB 



V X (v X B) - V X [77V X B 

4- vn{y X B) X B + ryA(V X B) 



(12) 



(e.g. |Wardle|2007| ). Others prefer to use separate equations 
for the charge carriers and the neutrals. In the latter 'multi- 
fluid' approach, when the ionised species are ions and elec- 
trons only (denoted by the subscripts 'i' and 'c', respec- 
tively), the induction eq uation becomes (e.g. |Konigl„1989l 



Balbus fc Terquem||2001 aano fc Stone|2002a| 



dB 



V X (v X B) - V X 



J X B _ [J X B)x B 



(13) 



Comparison of equations ( 12 1 and ( 13 1 yields the expres- 
sions given in ([6| to ([8| for the diffusivities. On the r.h.s. of 
equation ( 12 1, as well as in the second line of equation ( 13 1, 



the terms (from left to right) denote the inductive (I), Ohm 
(O), Hall (H) and ambipolar diffusion (A) contributions to 
the evolution of the magnetic field, respectively. 

The following parameters have typically been used to 
characterise the magnetic properties of the fiuid: 

(i) Field strength. This property is commonly measured 
either by the ratio of the Alfven speed to the isothermal 
sound speed (va/cs), or by the plasma beta parameter /3 = 
(2/7)03/11^, where 7 is the adiabatic index of the fluid. 

(ii) Field-matter coupling. The degree of coupling between 
the neutral matter and the magnetic field is, most generally, 
measured by a ratio of the type VL/D where V , L and D 
represent characteristic speed, length and diffusion scales of 
the flow, respectively. In MRI studies, it is appropriate to 
take V = va, L = wa/^ (the characteristic wavelength of 
MRI-unstable modes in ideal-MHD conditions) and D = ri±, 
the total perpendicular diffusivity [e.g. equation (C81]. The 
resulting ratio is the Elsasser number 



A; 



A. 



(14) 



This dimensionless number is formally deflned as the ratio of 
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Table 1. Typical dimensionless numbers used to characterise the degree of coupling between the magnetic field and the neutral gas: 
Magnetic Reynolds No., Lundquist No., and Elsasser No. For comparison, we also show the Reynolds No. (although this parameter 
combination does not involve the magnetic diffusivity and, therefore, does not measure the field-matter coupling). The fluid variables 
used to characterise the required scalings are the flow velocity (v), the Alfven velocity (da)i the kinematic viscosity {u) and the magnetic 
diffusivity rj. As usual, O is the Keplerian angular frequency and Z is a typical, unspecifled, length scale of the flow. 



Symbol 



Dimensionless 
number 



Ratio of 
physical quantities 



Scaling 
V L D 



Definition 
VL/D 



Re Reynolds 

Rem Magnetic Reynolds 

S Lundquist 

A Elsasser 



Inertial to viscous forces v 

Inertial to resistive time scales v 

Resistive diffusion to Alfven time scales da 

Lorentz to Coriolis forces v\ 



u vl/v 

r] vl/rj 

V vaI/v 
V± 



v%/rij_n 



the Lorentz force {J x B)/c to the Coriolis force oc p{v x CI), 
which results in an expression of the form v\/{LV)fl. Adopt- 
ing the magnetic diffusivity ri± as the typical magnitude of 



the LV factor in the denominator yields expression ( 14 1 



The Elsasser number is widely used in geophysics, as in a 
geodynamo process the magnetic field is thought to amplify 



until A becomes of order unity (see, e.g. Christensen 2010 
and references therein). This ratio was denoted x in the work 
of Wardle ( 1999 1, who expressed it as 



(15) 



where cr± = /(iivri^) and uic is the critical frequency 
above which flux-freezing conditions break down, so that 
the generic non-ideal MHD term in the induction equation 
dominates over the inductive term. When the field-matter 
coupling is characterised by the Elsasser number, values S> 1 
(<C 1) correspond to strong (weak) coupHng. Similarly, from 



equation ( 15 1 it is clear that when A ^ 1, the magnetic field 



is strongly coupled to the fiuid at frequencies of order Q, the 
frequencies of interest for the study of the MRI. 

Note that, in general, the Elsasser number is different from 
the combination vaI/t], which compares the resistive diffu- 
sion time Tpi oc to the Alfven time ta oc 1/va (e.g. the 
Lundquist number) and from the ratio vl/rj (the magnetic 
Reynolds number i?EM , where v is the fluid velocity) . How- 
ever, both terms have been used in the literature to refer 
to A. Table [l] summarises the definitions and typical scal- 
ings associated with the dimensionless numbers discussed 
above, namely, the magnetic Reynolds, Lundquist and El- 
sasser numbers. For clarity we also list the Reynolds num- 
ber, although this does not measure magnetic coupling, as 
it deals with the viscosity of the fiuid instead of its magnetic 
diffusivity. 

Finally, sometimes Cs is adopted as the characteristic ve- 
locity scale of the flow instead of ua, resulting in the ratio 
cl/r]Q. (e.! 



Fleming et al.[2000 1 



(iii) Diffusivity regime. This property characterises the 
importance of the different non-ideal MHD terms in the in- 
duction equation. The ratios of each of the diffusive terms to 
the inductive term are typically used, which in a multifluid 
formulation are expressed as (e.g. Balbus &: Terquem„200l| 
Sano fc Stone|2002a[ ): 



/ 



r]Q. 



(16) 



and 



H 

T 

A _ 

1 " 



_ X 

li Pi- 



rn 



(18) 



Note that the i nver se of the normalised Ohm resistivity 



f] in equation ( 16 1 has been referred to as the magnetic 



Reynolds number. However, as discussed above, this param- 
eter combination is akin to the Elsasser number, with the 



difference that in the definition of equation (141 we used the 



total perpendicular diffusivity t^j^, instead of r\. In equation 

(TtI, 



X 



cBn 

2-nenEv\ 



(19) 



is the so-called "Hall parameter", not to be confused with 
its namesake, the ratio of the gyrofrequency to the collision 
frequency of ionized species j (here either ions or electrons) 
with the neutrals, given by 

ft-^^. (20) 
rujC-yjp 

Similarly, the relative importance of the non-ideal MHD 



terms can be measured by the following ratios (e.g. Balbus 
fc Terquem|2001||Sano fc Stone|200"2al ) 

H_ _ eB__}_ ^ p _ ^ 

and 

A eB 



(21) 



Using the tensor diffusivity notation, we can write 
+ A_ ripQ 



H _ ?7H» 



(22) 



(23) 



(24) 



and the relative magnitudes of the non-ideal MHD terms 
are, simply 

H_ _ rjn 

O ^ V 

and 

A _ 77A 
H 77H ' 



(25) 



(26) 
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Dissipation 
suppresses 
MRI 



~ldeal-MHD 
MRI 




H/i = {m^)/vl 



Figure 9. Regions of parameter space where the MRI is expected 
to grow, and dominant diffusion mechanism, in a (0+A)/I versus 
H/I (or, equivalently, the normalized Pedersen versus HaU diffu- 
sivities) plane. Note that for the vertical field geometry considered 
here, the ambipolar and Ohm diffusivities can be combined into 
a generalized 'Ohm-type' diffusivity (see text). In the lower-left 
(ochre) panel, magnetic diffusion is weak and the instability grows 
at a rate comparable to the ideal-MHD rate (~ O). Conversely, 
in the upper-left (grey) and lower-right (purple) panels, only one 
diffusivity term (Ohm-type and Hall, respectively) dominate over 
the inductive term in the induction equation. Ohm-type diffusiv- 
ity suppresses the instability in the top-left quadrant (see foot- 
note|4] however), but since Hall diffusion can also be destabilising 
(Section [2]| the instability can still grow in the lower-right panel. 
Finally, both diffusivity components are important in the upper- 
right quadrant of the figure. In the top (grey) portion of this panel. 
Ohm-type diffusion overcomes Hall diffusion [H/(0-|-A) < 1] and 
the MRI is expected to be suppressed. Hall diffusivity, however, 
is dominant [H/(0-|-A) > 1] in the bottom (purple) portion of 
the panel. In this region of parameter space, H/(0-|-A) and H/I 
are both > 1 and the instability may proceed despite A being less 
than unity. 



From expressions ( 23 1 and ( 24 1 it is clear that the Elsasser 
number is a measure of the ratio of the inductive term to 
the total non-ideal MHD terms in the induction equation. 
Furthermore, for the vertical field geometry adopted here, 
the ambi polar diffusivity acts as a f ield-dependent resistiv- 



itjg(s 



Balbus & Terquem 20011. This property will be 



used in the next section to treat the Ohm and ambipolar dif- 
fusivities as an 'Ohm-like' term when comparing the impact 
of the 'Ohm-like' and Hall resistivities on the MRI. 



4.2 Criterion for Hall diffusion to affect the MRI 

We now use the ratios given in the previous subsection to 
constrain the region of parameter space where the MRI is 



^ In other words, for this field geometry the current density J is 
perpendicular to B and the term {JxB)xB = {JB)B-B^J = 



0.8 



0.6 



S 0.4 



0.2 





0.1 1.0 10 



2 4 
Log (H/0 = X/27]) 



Figure 10. Growth rate of the fastest-growing MRI-unstable 
mode as a function of the ratio of the Hall to Ohm diffusion 
terms in the induction equation [H/O = X/{2fj)]. Calculations 
correspond to a stratified disc with /3 = {2/f)c'^/v'j^ = 3200 (or 
va/cs = 0.02). Each curve corresponds to a different value of 17"^ 
(= I/O): 0.1 (blue), 1.0 (red) and 10 (green). A solid (dotted) 
line is used to plot solutions for which the ratio H/I > 1 (< 1). 
Note that solutions for which Hall diffusion is dominant lie in the 
solid-line portion of the curves (H/I > 1) and to the right of the 
vertical dashed line (H/O > 1). Superimposed (filled circles) are 
the linear growth rates corresponding to the solutions presented 
in Fig. 9 of Sano fc Stone| 1,2002b) for X = A. Note that the cal- 
culation for fi~^ = 0.1 and X = 4 satisfies H/O = 0.2 and is not 
in the Hall-dominated region of parameter space. As a result, the 
MRI is significantly damped by Ohm diffusion. On the contrary, 
the instability can grow at about the ideal-MHD rate when H/O 
> 1 even when fj~^ = 0.1. 



expected to grow, as well as to determine a criterion for Hall 
diffusion to substantially modify its properties. First, how- 
ever, we discuss the results of the work of |Sano fc Stone] 
(2002b) on this topic. These authors characterised the Ohm 
and Hall terms by their magnitudes relative to the inductive 



term in the induction equation (e.g. via eqs. 16 and 17 1. Fur- 



B J (see equation dlSfc). 



thermore, they considered that Hall diffusion was dominant 
when the ratio of the Hall to inductive terms was larger than 
unity (H/I > 1). Fig. 9 of |Sano fc Stone] (|2002b ) shows the 
saturation level of the Maxwell stress (normalised by the ini- 
tial gas pressure), as a function of the initial (subscript '0') 
of the inverse Ohm resistivity (?7o~^) - called the magnetic 
Reynolds number - for different values of the initial plasma 
beta and Hall parameters (/3o and Xo, respectively). 

Their results indicate that when r]o~^ > 1, the nor- 
malised saturated value of the Maxwell stress is of the order 
of 0.1 (with some scatter depending on the adopted /3o and 
Xo), and it is fairly independent of the actual value of the 
resistivity. The actual saturated magnitude of the stress is 
larger in the models with Xo — 4 with respect to the ones 
with Xo = (e.g. no Hall diffusivity) or Xo ~ —2 (negative 
Hall diffusivity, when the magnetic field is counter-aligned 
with the disc angular velocity vector). On the contrary, when 
the initial inverse resistivity ifo~^ is less than unity, the sat- 
uration level of the Maxwell stress decreases by 1-2 orders 
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of magnitude with respect to the ?7o~^ > 1 case. This trend 
of the Maxwell stress at saturation with 770"^ seems to be 
unaffected by the presence (or magnitude) of the Hall dif- 
fusivity. These results have been interpreted to show that 
the Hall diffusivity does not change the critical (maximum) 
fj required for the instability to grow (r^crit ~ 1), but it en- 
hances the saturated level of the Maxwell stress by a factor 
of a few. 



Note that the calculations of Sano & Stone ( 2002b I 



probe the region of parameter space where Ohm and Hall 
diffusivity terms dominate over the inductive term {I/O — 
f]^^ < 1 and H/I = Xo/2 > 1). However, |Xo| < 4 in all the 
presented calculations. As a result, for 770 > 2, the ratio H/0 
= Xo/(277o) is < 1 and the dominant diffusion mechanism is 
Ohm-type. For the calculations with Xq = 4 and fj^^ = 0.1, 
in particular, the ratio H/0 = 0.2, which implies that the 
Hall term is too weak to overcome the damping effect intro- 
duced by Ohm diffusion. It is not surprising, therefore, that 
the large drop in the saturated value of the Maxwell stress, 
with respect to that associated with the solutions satisfying 
tiq^ > 1, is not significantly modified by the presence of Hall 
diffusion. Naively, and as confirmed by the linear analysis, 
the Hall effect should substantially modify the growth rate 
of the magnetorotational instability provided that Hall dif- 
fusion dominates over the inductive term (i.e. |Xo| >2) and 
also dominates Ohm diffusion (i.e. Xo>2r;o)- According to 
these criteria, the Hall term would be strong enough to mod- 
ify the presented results if Xq > 20 (for the fi^^ = 0.1 case), 
a value likely to be well beyond what has been computation- 
ally feasible so far. 

The above considerations are summarised in Fig. [Oj 
which shows the regions of parameter space where the insta- 
bility is expected to grow - and, if so, the dominant diffusion 
mechanism - in a [(0+A)/I - H/I] or, equivalently, the nor- 
malised Pedersen and Hall diffusivities, plane. As discussed 
at the end of Section [4. 1| we use a generalised 'Ohm' term, 
which should be understood to mean 'Ohm + ambipolar'. 

Diffusivity effects are weak in the lower-left (ochre) 
quadrant of the figure, as the inductive term dominates over 
both the Ohm and Hall diffusion terms in this region of pa- 
rameter space. The MRI is then expected to grow at a rate 
not significantly reduced from the ideal-MHD rat^ Con- 
versely, in the upper-left (grey) and lower-right (purple) pan- 
els, only one diffusivity term is dominant over the inductive 
term (the Ohm and Hall term, respectively). When Ohm dif- 
fusion dominates, the instability's growth is damped, as this 
diffusion mechanism is always stabilising (See the discussion 
in Section |2] and Fig. [2|. However, the Hall term can be ei- 
ther stabilising or destabilising (see Figs.[3]and|4|, so when it 
dominates the instability can still potentially grow . Finally, 
the upper-right quadrant of the figure contains the region 
of parameter space where both the Ohm and Hall diffusiv- 
ity terms dominate over the inductive term. Although both 
diffusivity terms are important with respect to the ideal- 
MHD term in this region, only in the lower (purple) portion 

* Note, of course, that the boundaries between these regions are 
not sharp as this simplified diagram might suggest, and the MRI 
is expected to grow at reduced rate when rjpvj^/Q ^ 1 (e.g. near 
the border between the lower and upper- left panels of the figure) . 
This also applies along the diagonal line separating the Hall and 
Ohm-dominated regions in the upper-right panel. 



of the panel Hall diffusion is strong enough to potentially 
overcome the damping effect introduced by Ohm diffusivity, 
as only here the ratios H/(0-|-A) and H/I are both > 1. 
Note that this region of parameter space was not probed by 
'Sano & Stone (!2002b), as their solutions incorporating Hall 



diffusion lie in the upper (grey) portion of the panel, where 
Hall diffusion is weak in comparison to Ohm diffusion. The 
potentially destabilising effect of a sufficiently strong Hall 
diffusivity (the purple region of the upper-right quadrant 
in Fig. [9]l, therefore, remains to be explored by numerical 
simulations. 

In order to test these assertions, we computed the 
growth rate of the most unstable MRI mode in a strat- 
ified, geometrically thin and axisymmetric accretion disc, 



using the method described in Salmeron & Wardle (20031, 



to which we refer the reader for additional details. Our re- 



sults are depicted in Fig.[To]for P = 3200 (or va/c^ = 0.02), 
as a function of the ratio H/O (= X/2fj). Each curve cor- 
responds to a different value of the ratio I/O (= r;"^), as 
follows: fj-^ = 0.1 (blue), 1.0 (red), and 10 (green). The 
range of solutions shown with solid lines satisfy X > 2, so 
that the Hall term dominates over the inductive term (H/I 
> 1). The solutions depicted with a dotted line correspond 
to X < 2 (H/I < 1). Note also that to the left of the vertical 
dashed line the ratio H/O = X/2'i) < 1 and Ohm diffusion 
is the dominant non-ideal MHD term. 

From the considerations in the previous paragraphs, it 
is clear that the solutions that probe the Hall-dominated 
region of parameter space should lie in the solid segment of 
each curve and to the right of the vertical dashed line. Note, 
in particular, that even for fj~^ = 0.1, the instability grows 
at essentially the ideal-MHD rate (z-'max/fi ~ 0.75) if the 
Hall term is strong enough (H/0 > 1, or to the right of the 
vertical line). For comparison, we superimpose filled circles 
on the curves showing the linear growth rate of the most un- 
stable mode corresponding to X = 4 and fj~^ = 0.1, 1 and 
10 (the set of parameters tested by |Sano fc Stone]|2002b| . 
Note that the solution with X = 4 and ?7~ =0.1 lies to 
the left of the H/0 = 1 line, and Ohm diffusion suppresses 
the growth of the instability. We conclude, therefore, that 
the relative saturation levels of the MRI simulations pre- 



sented in Sano & Stone ( 2002b I are entirely consistent with 



expectations based on a linear analysis, and did not probe 
the regime H > {/, O, A} in which the hall term strongly 
modifies the linear MRI. 



5 APPLICATION TO PROTOPLANETARY 
DISCS 

We now illustrate the implications of Hall diffusion for the 
extent of MRI-driven activity in protoplanetary discs by ap- 
plying the linear analysis from appendix |B] to the minimum 
mass solar nebula at 1 au 

While a linear analysis cannot hope to give a good sense 
of the properties of the nonlinear turbulence the MRI drives, 
numerical simulations in the Ohm limit have shown that it 
does predict when such turbulence exists, and appears to be 
an excellent predictor of the vertical extent of dead zones 



in stratified simulations of protoplanetary discs (Turner & 
|Sano|2008l ). As noted in the previous section, the Hall-Ohm 
simulations of Sano & Stone ( 2002a|b I are also consistent 
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with the expectations based on the linear dispersion relation. 
Thus, we use the local dispersion relation to examine the 
role of the magnetic diflFusivity in determining the extent of 
MRI-active regions in protoplanetary discs. 

The tendency of the MRI to manifest on successively 
larger scales in the presence of increasing diffusivity is lim- 
ited by the finite thickness of the disc: the wave numbers of 
interest are bounded from below by requiring kh > 1 where 
h = Cs/Q. is the disc scale height. Thus we adopt as a local 
criterion for growth of the MRI that the dispersion relation 
(B6l yields modes with v > and kh > 1. 



From this we 

infer that the diffusiviti es m ust lie within the semicircular 

1 



locus obtained from eq (BIO I with u = Q and k = h 



I.e. 



3 cj^ 
4^ 



+ 



77pn 



l6 



1 + 



(27) 

where s = sign{Bz). This criterion applies to any near- Kep- 
lerian disc threaded by a vertical magnetic field. Its apparent 
simplicity belies the fact that the diffusivities are compli- 
cated (but calculable) functions of magnetic field strength, 
gas density, temperature, and the abundances of charged 
species. Note also that the ratio (?s/v\ plausibly ranges be- 
tween and the equipartition value 2 at the disc mid-plane 
but may be significantly higher away from the mid-plane if 
the field is anchored at lower heights within the disc. 

We specialise to protoplanetary discs by adopting the 
ionisation models of Wardle (20071 for a standard mini- 



mum mass solar nebula disc ionised by cosmic-rays and 
stellar x-rays at 1 an from the central star. These models 
assume a column density of 1700 gcm~^ and temperature 
T = 280 K independent of height, with x-ray ionisation rate 



computed by Igea & Glassgold ( 1999 1 and a standard in- 
terstellar cosmic-ray ionisation rate of 10~^^ s~^ that is 
exponentially attenuated with depth as exp(— E/96gcm~^). 
We characterise the grain population by assuming a single 
1 /xm radius and varying the dust-to-gas mass ratio, crudely 
mimicking the effect of the settling of grains to the disc mid- 
plane. 

The abundances of charged species for dust-to- gas mass 



ratios of 10 , 10 and are presented in Figure 11 The 



ionisation rate as a function of depth is plotted in the lower 
panel - interstellar cosmic rays are the dominant source be- 
low 2 scale heights, above this ionisation by stellar x-rays 
dominates. In the no-grain case, electrons and metal ions 
are the dominant charged species because the metal ions 
have the smallest recombination rate coefficient. The top 
panel shows the effect of adding a population of 1 /im-radius 
grains with total mass 1% of the gas mass: grains acquire a 
charge via sticking of electrons and ions from the gas phase. 
Above z/h « 2.5 the grain charge is determined by the 
competitive rates of sticking of ions and electrons, with the 
Coulomb repulsion of electrons by negatively-charged grains 
offsetting their greater thermal velocity compared to ions 
( |Spitzer|[794T] |Draine fc Sutin||1987[ ). This leads to a gaus- 
sian grain charge distribution with mean charge (in units 
of e) {Zg) ~ —AakT/e^ « —67 and standard deviation 
Ri {Zg)^^'^ « 8. Most recombinations still occur in the gas 
phase. The abundance of metal ions and electrons is reduced 
over the grain-free case for < 5 where the charge stored 
on grains becomes comparable to the electron abundance. 
For z/h « 2-2.5, the abundances of ions and electrons have 
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Figure 11. Fractional abundances of charged particles in the 
minimum mass solar nebula at 1 au from the Sun as a function 
of height above the mid-plane, assuming that grains have radius 
1 ^m. The upper, centre, and lower panels correspond to dust to 
gas ratios 10~^, 10~*, and by mass, respectively. The black 
curve in the lower panel shows the height-dependence of the ion- 
isation rate per hydrogen nucleus due to stellar x-rays and inter- 
stellar cosmic rays (see text) assumed for all three models. The 
other curves give the fractional abundances of electrons (red), 
light ions (H+, Hg", He+, 0+) representative molecular (m+), 
and metal (M+) ions (blue), and grains {green, labelled by charge 
state). 
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declined to the point that the majority of electrons stick to 
grain surfaces before they can recombine in the gas phase, 
and most neutralisations occur when ions stick to negatively 
charged grains ( Nakano fc Umebayashi|1980 1. Closer to the 
mid-plane, the ionisation fraction is so low that most grains 
are lightly charged, so that Coulomb attraction or repulsion 
of ions and electrons by grains is negligible. Then ions and 
electrons stick to any grain that they encounter, with re- 
combinations occurring on grain surfaces. The middle panel 
of Fig. [ll] shows what happens if 99% of the grains are re- 
moved (e.g. by settling to the mid-plane). The capacity of 
the grain population to soak up electrons from the gas phase 
is reduced a hundredfold, and the height below which grains 
substantially reduce the free electron density below the ion 
density moves downwards to the lowest scale height. 

Unlike Ohm resistivity, the Hall and ambipolar diffu- 
sivities depend on the magnetic field strength as well as the 
charged particle abundances, so we consider field strengths 
ranging between 10~'^-10'^ G, encompassing the plausible 



'dust l^gas = ^ 



range of values in the solar nebula (see Wardle 20071. For 
each choice of magnetic field we use the ionisation calcu- 
lations to compute the diffusivities as a function of height. 
Typically, ambipolar diffusion dominates at the surface and, 
for weak magnetic fields. Ohm diffusion dominates near the 
mid-plane. Between these regimes Hall diffusion dominates 
( Wardle|2007 l. With the diffusivity profiles in hand, we ap- 
ply eq (27 1 to identify the fastest growing MRI mode that 



has fc/i > 1 at a given height; the results f or d ust -to- gas 
are displayed in Figs 



12 



and 



13 



mass ratios and 10 
spectively. 

First consider the case in which grains are absent 
(Fig. |12[ ). The bottom panel shows the fastest growth rate 
(subject to kh > 1) as a function of magnetic field strength 
and height if Hall and ambipolar diffusion are neglected and 
only Ohm diffusion is included. In this case the departure 
from ideal MHD is measured by rjQ,/v\ (the inverse of the 



Elsasser number; see section 4.11, which strongly decreases 



with height as the fractional ionisation and Alfven speed 
both increase sharply away fro m t he midplane. For the ion- 
isation profile plotted in Fig. 11 it turns out that riQ,/v\ 



is small near the surface (so that ideal MHD holds) and 
large near the mid-plane (so that Ohm damping is severe). 
As a result, near the surface the largest achievable growth 
rate is close to the ideal value 0.75 Q., but declines rapidly 
below the height where r}Q/v\ ~ 1. The height of the tran- 
sition between these regimes declines with increasing field 
strength, simply because 7)0. /v\ cc . A second consid- 
eration is that the range of MRI-unstable wave numbers is 



bounded above by fee (see eq B12 I, and this must be larger 
than l//i if any unstable modes are to exist with fc/i > 1. It 
is this criterion that provides the upper and lower envelopes 
to the unstable region in Figs. [12] and |13| Near the surface 
where r]Q,/v\ ^ 1, fee is approximately y/ZO/vA, and as 
va rapidly increases with height the unstable range shifts to 
wave numbers with kh < 1 and it is no longer possible to find 
unstable modes that fit within a scale height. This occurs 
when the magnetic pressure roughly exceeds the gas pres- 
sure. By contrast, close to the mid-plane where tiQ/va ^ 1, 
fee ~ VSva/v and the wave numbers are pushed out of the 
relevant range by the increasing diffusivity and the declining 
Alfven speed as the mid-plane is approached. 

Next, we add the remaining magnetic diffusion terms. 
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Figure 12. Colour shading shows the growth rate of the fastest 
growing local MRI mode at height z above the mid-plane in the 
minimum mass solar nebula at 1 au as a function of the strength 
of an initially vertical magnetic field B. Dust grains are assumed 
to have settled to the mid-plane (see lower panel of Fig. 
The vertical wave number k of the modes are required to satisfy 
kh > 1 where h is the disc scale height. The unshaded regions 
show stable combinations of height and field strength. The top 
and middle panels correspond to the cases where the initial field 
and rotation axis are parallel {B^ > 0) or antiparallel {B^ < 
0), respectively. The lower panel shows the effect of artificially 
suppressing ambipolar and Hall diffusion, so including only Ohm 
diffusion. In this case there is no dependence on the sign of Bz . 
The dashed line in each panel indicates the value of the local 
equipartition magnetic field as a function of height above the 
mid-plane. 
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ie. ambipolar and Hall diffusion. The latter depends on the 
sign of Bz, which we take to be positive or negative in the 
upper and middle panels of Fig. [12] respectively. Ambipo- 
lar diffusion dominates in the surface layers (see Fig. 5 of 
Wardle|[2007 1, and it acts just like Ohm diffusion when, as 
assumed here, the initial magnetic field is vertical tend- 
ing to damp the growth of the MRI and pushing the un- 
stable modes to longer wavelengths. Hall diffusion tends to 
dominate the lower layers, except for weak fields when Ohm 
diffusion dominates closer to the mid-plane. 

For Bz > Hall diffusion is destabilising and pushes 
the region where the fastest growth rate is close to 0.75 Q 
to greater depths. The rapid decline to low growth rates 



Mdust/Mgas = 10 2 



(ie. from green to blue shading in Fig. 12 I traces the transi- 
tion from Hall-dominated to Ohm-dominated diffusion. Be- 
low this. Hall diffusion acts to significantly extend the region 
of slow growth by modifying kc, and this extends all the way 
to the mid-plane for fields in excess of 10 mG. 

On the other hand, when Bz < (ie. B is antiparallel 
to the rotation axis) Hall diffusion tends to stabilise the 
MRI while extending the unstable wavelengths to shorter 
wavelengths (see Fig. [6|. As a result, super-equipartition 
fields are unstable near the mid-plane. We emphasise that 
irrespective of the magnitudes of the other diffusivities Hall 
diffusion is completely stabilising when rja > 2wx/f2.This is 
responsible for the sharp cutoff at the lower boundary in the 
middle panel of Fig. |12| 

When a full complement of 1 fim dust grains are present 
(i.e. dust-to-gas mass ratio 10~^; charged species as in lower 
panel of Fig. Ill, the diffusivities are greatly increased near 



the mid-plane because electrons are locked up by grains and 
rendered immobile. Fig. [13] displays the same trends as in 
the zero-grain case, but the MRI-unstable region is now re- 
stricted to the upper layers of the disc, and in all cases the 
bulk of the disc is stable to the MRI. While the differences 
between the three panels might appear less severe in this 
case, the strong density stratification means that there are 
orders of magnitude differences in the column density of the 
MRI-unstable region. 

Fig. [14] shows the magnetically-active column density 
as a function of field strength for different assumptions re- 
garding the diffusivity when 99% of grains are assumed to 
have settled (i.e. dust to gas ratio 10~*). The active column 
density varies by 1-2 orders of magnitude depending on the 
accuracy of the treatment of magnetic diffusion. First, ne- 
glect Hall diffusion and consider either Ohm diffusion alone 
or Ohm and ambipolar diffusion (i.e. Pedersen diffusion) op- 
erating in concert (red long-dashed curves). Ohm diffusion 
dominates ambipolar diffusion except for strong magnetic 
fields and low densities, so there is little difference between 
these two cases except for magnetic fields in excess of 1 G 
where the additional damping due to ambipolar diffusivity 
suppresses the MRI. Hall diffusion either increases or de- 
creases the active column density by an order of magnitude 
depending on whether the initial magnetic field is pointing 



^ This is not the case in more general geometries, e.g. ambipolar 
diffusion may be destabilising when the field has vertical and 
toroidal components and the wave number has radial and vertical 
components (see Kunz &: Balbus|2004 and Desch|2004 1 . 




N 
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Figure 13. As for Fig. |12| but now including a population of 
1 ^m radius grains with total mass 1% of the gas mass (see upper 
panel of Fig.[IT]l. 



up or down; the maximum unstable field strength also varies 
considerably. 

Fig. [15] shows how the active column depends on grain 
abundance and magnetic field orientation. In the absence 
of grains the entire disc cross-section is magnetically active 
for an upwardly-directed field in the range 0.01-5 G. On the 
other hand, when the field is pointing downwards, a minus- 
cule fraction of the disc is active for weak fields, but the 
active column increases rapidly and encompasses the entire 
cross section of the disc for field strengths in the range 20- 
80 G. As one would expect, grains sharply reduce the active 
column density because of the reduction in mobile charge 
carriers. The continued extreme sensitivity to field align- 
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B(G) 



Figure 14. Column density of the MRI-active region at 1 au in 
the minimum mass solar nebula as a function of magnetic field 
strength for different assumptions about magnetic diffusion (see 
text). Grains are assumed to have radius l/xm with 99% having 
settled to the disc mid plane (i.e. total mass only 10~^ relative to 
gas). Solid black and dashed blue curves are for the magnetic field 
oriented parallel or antiparallel to the disc's rotation axis respec- 
tively. Long-dashed red curves indicate the effect of accounting 
only for Ohm or Ohm+ambipolar (i.e. Pedersen) diffusion. 




B(G) 

Figure 15. As for Fig. |14| but for dust-to-gas mass ratio vary- 
ing from to 10"'^. Note that active column density is capped 
at the total column of the minimum mass solar nebula at 1 au, 
i.e. 1700gcm-2. 



ment demonstrates that Hall diffusion still plays a critical 
role in determining the extent of magnetic activity regard- 
less of the grain abundance. 

The reduction in the active column when Bz < occurs 
because Hall diffusion stabilises the disc against the MRI 
when srjH < —2v\/Q. When the grains and ions are strongly 
coupled to the neutrals and the electrons are coupled to the 
magnetic field, the Hall diffusivity is given by eq ([8| and 
we obtain a field-dependent criterion for stability on the 
fractional ionisation 



^ < 
n^r ~ 2 



1 lAruHcn 1.5 X 10" 



(M/Me 



\l/2 



eB 



B{G) 



r/au)^/^ 



(28) 



boundary of the magnetically active region when B^ < 0. 
Note that once the grain abundance is low enough not to 
affect the ionisation fraction, the lower boundary reaches 
the midplane and the active column density saturates. 

These results suggest that Ohm estimates of the column 
density of the magnetically-active layers in protoplanetary 
discs are in error by about an order of magnitude, system- 
atically under- or overestimating the active column if the 
magnetic field is directed upwards or downwards, respec- 
tively. 

This conclusion rest on several assumptions. Our ne- 
glect of stratification in the linear analysis is unlikely to be 
serious given that in the Ohm limit this approach yields 
excellent predictions of the extent of the depth of the turbu- 
lent layers ( Turner fc Sano|2008 1. The restriction to a simple 
and somewhat degenerate geometry - i.e to vertical fields 
and wave vectors - is perhaps more suspect. However, when 
ambipolar diffusion is unimportant, a toroidal component of 
the magnetic field and a radial component of the wave vector 
can easily be accommodated: a rescaling of the full disper- 
sion relation yields the identical dispersion relation (Pandey 
& Wardle, in preparation). The situation is less clear when 
ambipolar diffusion is involved, as in this case it stabilises or 
destabilises the MRI | |Kunz Balbus||2004[ |Desch||2004[ ) in 
an analogous manner to the Hall contribution for the verti- 
cal field considered here (Pandey & Wardle, in preparation). 
However, at ~ 1 au, ambipolar diffusion only dominates near 
the disc surface where the ionisation fraction is so high that 
diffusion is small in any case. Finally, although the grain 
model is crude, the key property of the grain population is 
its capacity to soak up electrons, which is proportional to 
/ an{a)da ( | Wardle|200"7| . This can be used to scale our re- 
sults to any grain size distribution. The real uncertainty is 
the small-radius tail of the grain size distribution which may 
still be sufficient to reduce the electron and ion fractions at a 
few scale heights despite the tendency of grains to aggregate 
and settle towards the mid-plane. 

The robust point is that the relative magnitudes of the 
diffusivities are independent of the ionisation fraction, being 
controlled by the ratio B/nn, and that Hall diffusion dom- 
inates the other mechanisms over much of protoplanetary 
discs |Wardle fc Ngi|1999 : Balbus fc Terquem|[2001 Sano"fcl 



Stone 



2002a I. As X-rays and cosmic rays both penetrate to 



This critical fractional ionisation is typically larger than the 
corresponding criterion for Ohm damping, so sets the lower 



the Hall-dominated region just below the disc surface, it is 
absolutely essential that Hall diffusion be included in mod- 
elling the extent of MRI-driven turbulence in protoplanetary 
discs. The real uncertainty is that, notwithstanding the pi- 
oneering effort of Sano fc Stone ( 2002a|b l, the saturation of 
MRI-driven turbulence in the Hall-dominated regime as yet 
remains unexplored. 



6 SUMMARY & CONCLUSIONS 

In this paper we re-examined the role of Hall diffusion in 
suppressing or enhancing MRI-driven magnetic turbulence 
in Keplerian discs. 

We first undertook a local, linear analysis of the magne- 
torotational instability, for simplicity restricting our atten- 
tion to a vertical, weak magnetic field subject to axisymmet- 
ric perturbations with a purely vertical wave vector. While 
this is not new, our approach and presentation differed from 



© 2011 RAS, MNRAS OOO.ppOl 



16 Mark Wardle & Raquel Salmeron 



previous analyses ( Wardle|1999 Balbus fc Terquem|2001[ ) in 
two critical ways. First, we characterised the departure from 
ideal MHD using magnetic difTusivities rather than conduc- 
tivities or characteristic frequencies, allowing us to make a 
clearer connection with the E X B field-line drift implicit 
in the diffusive MHD induction equation. Second, we pre- 
sented a clear perspective of the dependence on the Hall and 
Pedersen (Ohm-|-ambipolar) diffusivities by examining how 
the properties of the instability vary over a Hall-Pedersen 
diffusivity plane (see Fig. [5|. This clearly delineates when 
the range of unstable wave numbers is finite or infinite, and 
whether the fastest growth occurs for finite wave number 
or asymptotically as fc — >■ oo. We also discussed the limit- 
ing forms of the dispersion relation, making a connection 
with the diffusive plane-parallel shear instabilities of |Kunz| 
(200811. 



Next, we reviewed the alternative parameterisations of 
non-ideal MHD that have appeared in the literature, and 
emphasised that existing simulations of the nonlinear devel- 
opment and saturation of the instability in the Hall-Ohm 
case ( |Sano fc Stone]|2002a|b"l ) are consistent with expecta- 
tions based on the simple linear analysis and have not yet 
probed the Hall-dominated regime characteristic of proto- 
planetary discs. 

Finally, we illustrated the critical effect of Hall diffusion 
on the size of dead zones in protoplanetary discs by applying 
a local criterion for growth of the magnetorotational insta- 
bility to a simple model of minimum-mass solar nebula at 
1 au , including x-ray and cosmic-ray ionisation and a pop- 
ulation of 1 /im grains. 

Our key results can be summarised as follows. 

(i) The radial diffusion of perturbed field lines through 
the fiuid is directly related to the criterion for marginal sta- 
bility. The MRI is suppressed if the radial drift of field lines 
against the infall of the fiuid is sufficient to restore the field 
to its equilibrium radial position. 

(ii) The behaviour of the MRI for our adopted geometry 
is determined by the dimensionless Pedersen and Hall dif- 
fusivities ripfl/vA and srin^/v^, where s — sign{Bz). For 
srju < —2v\/Q. Hall diffusion suppresses the MRI irrespec- 
tive of the value of jyp. For sijn > —2v\/Q. the rjp — STjn half 
plane is divided into three regions in which either (I) there is 
a finite range of unstable wave numbers with a fastest grow- 
ing mode; (II) instability for all wave numbers, with a unique 
fastest growing mode and a slower asymptotic growth rate as 
fc — >■ oo, or (III) instability at all wave numbers with fastest 
growth asymptotically achieved as fc — > oo (see Fig. [sjl . 

(iii) For fixed ryp, the maximum growth rate increases 
with increasing S77H, from at sr/nf^/wA = — 2 to the maxi- 
mum 0.75 Q, as s?7h — >■ +00. For fixed 77H with st^h > —2v\Q, 
the growth rate is a maximum for j^p = and declines as 
77P — >• 00 (see Fig. [6|. 

(iv) In the highly diffusive limit the instability reduces to 
the Hall-diffusion-driven instability in plane-parallel shear 



fiow discussed by Kunz (20081. Diffusion is so severe in this 



limit that the perturbations in the fiuid velocity do not affect 
the field evolution, which is driven purely by diffusion and 
Keplerian shear. Our restriction to vertical initial fields and 
perturbation wave numbers enabled us to extend the results 
of |Kunz| ( |2008[ ) to include the damping by Pedersen diffusion 
and to derive pleasant analytic expressions for the upper 



cutoff to the unstable wavenumber range and the growth 
rate and wave number of the most unstable mode (see eqs 
C7| C9 and CIO respectively, and also Fig. [8|. 

(v) We argued that simulations of MRI-driven MHD tur- 
bulence in the presence of Hall and Ohm diffusion ( Sano & 
Stone|2002a|b \ have not yet probed the "deep" Hall regime 
*'?H ^Vp ^ ^a/^! where the linear analysis suggests that 
Hall diffusion allows the instability to proceed when it oth- 
erwise would not. 

(vi) We found that at 1 au in the minimum-mass solar 
nebula. Hall diffusion changes the magnetically active col- 
umn density by an order of magnitude. This change is either 
an increase or decrease depending on whether B is parallel or 
antiparallel to the rotation axis, respectively. Hall diffusion 
likely plays a critical role in determining the radial extent of 
dead zones and the thickness of magnetically active layers 
in protoplanetary discs, and estimates based on damping by 
Ohm diffusion are probably wildly inaccurate. 

The simplifications adopted in our analysis engender three 
significant uncertainties in our conclusions that are worth 
some final discussion. 

First, the restriction to the stability of vertical mag- 
netic fields to perturbations with vertical wave vectors does 
not capture the destabilisation by ambipolar diffusion that 
arises when toroidal field and radial wave vector components 
are also present (IKunz fc Balbus||2004 Dcsch 2004). How- 



ever, we have captured the analogous Hall diffusion-driven 
destabilisation of the MRI, and as Hall diffusion typically 
dominates this restriction is unlikely to have a great impact. 

Second, while the use of the linear analysis to predict 
the boundary of the manifestly nonlinear active region, ap- 



pears to be justified for the Ohm case ( Turner fc Sano 2008 1 , 



it is not known whether this applies in the Hall-dominated 
regime that we tout here. 

Finally, we caution that the MRI may simply be irrel- 
evant in protoplanetary discs. A minor population of small 
dust grains would remove so many electrons from the gas 
phase that the MRI-active column density becomes so small 
as to be irrelevant. Instead, magnetic activity - if any - may 
be due to fields lying within the disc and varying on length 



scales of order r ( Turner fc Sano||2008 1 or a strong poloidal 



magnetic field brought in during the formation of the disc. 
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APPENDIX A: FORMULATION 

Here we give the MHD equations describing a non-self- 
gravitating disc orbiting a point mass M. We adopt cylin- 
drical coordinates (r, 0, z) centred on M with the disc mid- 
plane corresponding to 2: = and the disc angular velocity 
vector parallel to the z-axis. We write the fluid equations in 
the inertial "laboratory frame" with the fluid velocity writ- 
ten as D -I- i;k, so that v denotes the departure of the flow 
from the Keplerian velocity field 



GM - 

9, 

r 



(Al) 



which we eliminate in favour of the Keplerian angular ve- 
locity 



Vk . 
— z 
r 



(A2) 



using the identities V ■ uk = and V X ~ This 
approach avoids the complexity that would be introduced 
by transforming to a local non-inertial frame or adopting a 
shearing-sheet approximation, neither of which are needed 
in our simple linear analysis in the next section. 
The continuity equation becomes 



(A3) 



The advective term containing corresponds to {vk ■ V)p 
and represents the azimuthal advection associated with Kep- 
lerian rotation - analogous terms will arise in the momentum 
and induction equations. The momentum equation. 



2Qv^r - 



V$- -VP + 



J X B 



(A4) 



P cp 

also picks up f2 -bearing terms that account for centripetal 
acceleration and the angular momentum loss implicit in ra- 
dial motion if the azimuthal speed is Keplerian. The gravi- 
tational potential is 

GM 



(A5) 



\^ r^ + z'^ 

and we adopt an isothermal equation of state: 

P = pc! , (A6) 

where Cs is the isothermal sound speecj^ The current density 
satisfies Amperes Law: 



J = —V X B, 

An 

and the magnetic field is, of course, solenoidal, 
V- B = 0, 



(A7) 



(A8) 



and evolves according to the induction equation (A9L The 



qualitative discussion of the effect of field diffusion and the 
magnetorotational instability in i|2] suggests that it will be 
useful to recast the induction equation in a form that makes 
explicit the drift of the magnetic field through the fluid, i.e. 



vb) X B -?7(V X B) 



(A9) 



^ This will not really be needed as the mode of interest does not 
involve pressure fluctuations. 
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where 



i;b = r]p 



(V X ■B)x X B 
B 



(V X B) 



B 



(AlO) 



Note that when vanishes, ambipolar and Ohm diffusion 
behave identically, appearing only together in sum as 77P. 
This is the case in the linear analysis of the MRI presented 
here. 

Note that the local dissipation rate associated with 
magnetic diffusion is 



J-E' [rjjfi+vpjl 



(All) 



and so, as is well known, Hall diffusion has no associated 
dissipation. When vanishes, ambipolar and Ohm diffusion 
behave identically, appearing only together in sum as rjp ; in 
particular this is the case in the linear analysis of the MRI 
presented here. 

Despite the appearance of 57-bearing terms, the equa- 
tions we have derived apply to an arbitrary fluid flow - we 
have not yet assumed that the flow is close to Keplerian, 
or even disc-like, but have simply written the fluid velocity 
in the MHD equations as v + vk- Of course, this form of 
the equations is only really useful for nearly Keplerian flows 
with <^ Vk- In the next section we shall specialise to a 
near-Keplerian disc equilibrium state with a vertical mag- 
netic field, and consider perturbations with a vertical wave 
vector. 



APPENDIX B: LINEAR ANALYSIS 

We consider a small region of an axisymmetric, geometri- 
cally thin and nearly Keplerian disc, threaded by a vertical 
magnetic field, with sound and Alfven speeds (c^ and va) 
at the mid-plane that are both small compared to the local 
Keplerian speed vk- We assume that radial gradients are on 
the scale of r and neglect vertical stratification of the initial 
equilibrium state, so that our analysis only holds near the 
mid-plane, at heights z <C C s/0.. Then, we may neglect the 



term rO, r — V"I> in eq (A4l, and the remaining radial and 



azimuthal derivatives appearing in eqs (A3|-(A9|, leaving 
only partial derivatives in t and z. 

The initial state has 1; = (i.e. is in Keplerian rotation) 
with a uniform density, pressure, and vertical magnetic field 
B = sBz (where s — ±1). We linearise the equations around 
this state and seek solutions for axisymmetric perturbations 
of the form exp(i^t — ikz). The equations for perturbations 
in density, pressure and Vz form a separate system that de- 
scribes vertically-propagating sound waves. The system of 
linear equations in the remaining perturbations involve fiuc- 
tuations in the r and (j> components oi B,v and ub. We treat 
these as two-component vectors in the equations below. 

In the following equations, we recast physical quantities 
in dimensionless form by adopting Q. and va as the units of 
frequency and velocity, respectively. For the remainder of 
this section, t^p, sryn, v and k denote r]pQ,/vA, rjn^/vA, v 
and kvA/^, respectively. Then the linearised momentum, 
field-line drift, and induction equations are 



Sv —ik 

" (l + u'^) 



SB 

B 



(Bl) 



VA 



-ik 



riP srju 
-srin ijp 



SB 



and 



B 



VA 



VA 



(B2) 



(B3) 



Using eqs (Bl I and (B2 I to substitute for Sv and Svb in eq. 

iSB, 

(B4) 



B3[ ) yields the relationship between the components of SB, 
sriH + 2/(1 + i^^^ 



5Br 



rfp + vA 



SB^ 



where 



1 



and the dispersion relation 

ak* + bk'^ + c^0 

where 

a = (1 + i'^){ri1> + r]u) + ^sr]H + 2iyrjp + 1 , 



(1 + u^) {2vriP 



+ 2u 



(B5) 

(B6) 

(B7) 

(B8) 

(B9) 
( Wardle|[T999| n 

While this form of the dispersion relation 
is convenient for determining the run of u with k for a given 
choice of the diffusivity components, it does not provide an 
overview of the dependence of the characteristic properties 
of the instability - the range of unstable wave numbers, 
maximum growth rate and corresponding wave number - 
on magnetic diffusion. To this end, we recast the dispersion 
relation into a form that emphasises its dependence on the 
diffusivities 77H and 77P: 



1 + 



3A 
4 



+ {rfp + vA) = ' 



(BIO) 



Conveniently, a given growth rate and wave number corre- 
sponds to a circular locus in the S77H — rjp plane, centred 
at (3A/4 - 2/(1 + i/^), -vA) with radius 3A/4, making this 
form amenable to graphical analysis. 

Using this expression it is straightforward to show that 
the maximum possible growth rate is = 3/4 (in units of 
Q,) and that this occurs only when rjp = 0. To see this, first 
not e tha t 77P > 0, and A> Q. Then, for u the LHS of 



eq (BIO I is > u A , with equality attained only if ?7p = 



and sr/H = SA/A — 2/( 1 + 1/ ). This statement must also hold 



for the RHS of (|B10}, so we conclude that i/^ < 9/16 with 
equality only holding when 77P = and srjn = 3/4fc^ — 5/4. 

Consider now the dependence of the instability on the 
Pedersen and Hall diffusivities rjp and rju (see Fig.jsJ. First, 
note that the ideal-MHD limit holds at the origin (i. e. rju ~ 
rjp — 0), the Hall MHD limit holds along the horizontal 
axis [rjp = 0) and the Ohm or ambipolar diffusion limits 
hold along the vertical line (rjn — 0). Recall also that only 
the half plane r;p > is physically relevant. Inspection of 
the dispersion relation, eq. (BIO I, shows that all modes are 
stable for srju < —2, and in the unstable region {srju > —2) 
there is, at most, a single unstable mode for a given choice of 



^ There is a typographical error in eq (14) of |Wardle| ( |1999| | — the 
RHS should be preceded by a minus sign. 
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wave number k. The run of growth rate with wave number, 
i.e. ^{k), for a particular choice of rjn and rjp, can be found 
by directly solving (B6l for v. In practice, it i s easier to 
choose u and solve the quadratic equation |B6l for k^. It 



turns out that there are three distinct forms of the resulting 
fik) curve, corresponding to regions I, II and III in the S77H- 
rjp plane, as illustrated in Fig. [5] 

In region I, which lies outside the semi-circular locus 



vl + {srni + 5/4)^ = 9/16 , 



(Bll) 



the range of unstable wave numbers extends from fc = up 
to a maximum value 



kc — 



srjH + 2 



2 »7p + {sva + 5/4)2 _ g/ig 



1/2 



(B12) 



found by setting v = and k = kc in the dispersion relation 



(BlOl. The form of ^{k) in this region is illustrated in inset 
I of Fig. |5] The maximum growth rate up c orresponds to 
a repeated root for k^ in the quadratic ( B6 1 , so it can be 



found by setting the discriminant to zero. This yields 



S77H = 



9 



VP 



1 + ^0' 



(B13) 



so contours of constant i/q in region I are straight lines. The 
corresponding wave number, ko, is then given by the ratio 
—2c/b, or alternatively — ti/2a, obtained from eqs (|B7|-([B9| 
with V set to vq, i-e. 



ko 



(1 + ui) (2i'o77P 



2^2 



1/2 



(B14) 



Region I encompasses the limits of ideal-MHD, the Ohm 
and/or ambipolar diffusion limit, as well as the Hall limit 
(i.e. -qp = 0) for srjn > —0.5. Elsewhere, i.e. within the 
semicircle bounded by (Bll 1 and the sr/n-axis {rjp = 0), all 
wave numbers are unstable. The interior of the semicircle 
is further subdivided into regions II and III, depending on 
whether the maximum growth rate is attained at a finite 
wave number or asymptotically as fc — > 00, as illustrated in 
insets II and III in Fig. [5] respectively. 

In region II, uq and ko still satisfy eqs ( |B13[ ) and ( |B14[ ) 
just as in Region I. The inner boundary of region II occurs 



where fco just becomes infinite, ie. i/q satisfies eq (B13I and 



simultaneously the denominator of eq (B14I is zero. These 



two conditions yield a parametric solution for the locus sep- 
arating regions II and III, 



S77H 



VP = 



2(9 + iui) 



(l + i.,f)(9+ 16^/2) 



^/o(9- IGi^g) 
(l + ;/2)(9 + 16i/o') 



(B15) 



which traces out an arc from (ryp, stjh) ~ (0, —2) to (0, — |) 
as vo runs from to 3/4 (the blue locus in Fig. [sjl. 

In region III, all wave numbers are unstable, with the 
fastest growth occuring in the limit fc — > 00, yielding from 



(BlOl 



sr^H + 



5/4 
1+^ 



+ Up + 



1^0 



3/4 
1 + ^0' 



(B16) 



Thus in region III, contours of constant fo trace out seg- 
ments of non-concentric circles running between the sr^H-axis 
and the boundary with region II. 

Having delineated these three regions, we now consider 
how the critical wave number k^, fastest growth rate vq and 
corresponding wave number ko vary across the entire 77p-S7;h 
plane. Contours of constant kc axe semicircles, as plotted in 
Fig. [7] While the range of unstable wave numbers is reduced 
for large values of S77H and r/p , as one might expect, the range 
is not maximised in the ideal limit (ie. at the origin) but in 
regions II and III, bounded by the fcc = 00 contour. 

Turning now to the fastest growing modes, the maxi- 
mum growth rate is given either by (B13l in regions I and 



II or by (BI61 in region III, and the corresponding contours 



are plotted in Fig. [6] The growth rate increases clockwise, 
from Ofl along the vertical line srju ~ —2 up to 0.75 f2 for 
the horizontal line r;p = for srju > —4/5. In the absence 
of Hall diffusion, the maximum growth rate i/q declines with 
increasing (Ohm and/or ambipolar) diffusivity (e.g. mov- 
ing vertically upwards through the srju ~ point in the 
horizontal axis), with uo ~ '^Vp^ for 3> 1. The most im- 
portant effect of Hall diffusion, apparent from Fig. [6] is that 
the growth rate of the MRI exceeds 0.3Q for srjn ^ rjPi even 
for arbitrarily large rjp. More generally, the addition of Hall 
diffusion at fixed rjp increases the growth rate if S77H > 
and d ecreases it when srju < 0. For large values of t^p, eq 
(B13I shows that STju/riP ~ 24t'o/(9 — IGuq). It is this fact 
that has the potential to modify the extent of dead zones in 
protoplanetary discs, as we explore later in !|5] 

The wave number of the fastest growing mode (blue 
contours in Fig. [6| decreases as the diffusivity is increased. 
Again, the contours are not arranged so that the highest 
wave numbers occur in the ideal-MHD limit, but to the 
sriH < side, within the boundary between regions II and 
HI (traced by the ko — 00 contour). 

Overall, these patterns place the ideal. Ohm (or am- 
bipolar) and Hall regimes in context, and for the first time 
we see an overview of the effect of magnetic diffusivity on the 
linear MRI. In particular, there is nothing special about the 
Ohm/ambipolar limit, e.g. the behaviour of the instability 
in the presence of diffusion is not qualitatively different for 
s»7h = 2 vs t;h = 0. Even the ideal-MHD limit does not stand 
apart as remarkable, although it still holds a special place 
conceptually because flux freezing holds and it is easier to 
think about. What does stand out is the part of the plane in 
the lower left, regions II and HI, characterised by high wave 
numbers and the spraying out of the growth contours. 

Overall we note that increasing rjp decreases the max- 
imum growth rate and the characteristic wave numbers, 
whereas increasing srjn above —2 increases the maximum 
growth rate and may either increase (when srjH + 1 ^ »7p) or 
decrease (when srjH + 1 ^ '7p) the corresponding wave num- 
ber. 



APPENDIX C: LIMITING CASES AND 
DIFFUSIVE INSTABILITIES 

In this section we consider the interesting limiting cases of 
our analysis. As a preliminary, we substitute equations (|B1[) 



and ( B2 \ into ( B3 1 to express the linearised induction equa- 
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tion in the form 



V I 

\ - 

VP sVa 

-S77H VP 



SB ^ 0. (CI) 



The three terms in this expression represent the effects of 
Keplerian shear, fluid displacement, and magnetic diffusion 
on the magnetic held perturbations. Each of three Hmiting 
cases can be obtained by neglecting one of these terms. To 
obtain the criteria for each limit, we note that i> < \, and 
so the three terms have order of magnitude 1, , and k'^ri±, 
respectively. The three cases of interest are as follows. 



Ideal MHD When k^ri± < fc^ ~ 1, i.e. r]± < 1, fe^ - 1, 
the third term is negligible, and held evolution is determined 
by shearing of the held and the response of the fluid to 
magnetic stresses. This limit applies in the neighbourhood 
of the origin in Fi g.js) an d (of course) recovers the results of 
[Balbus fc Hawley] ( |1991[ ) limited to a vertical fleld and wave 
vector and neglecting buoyancy. 



Cyclotron limit The flrst term is negligible when 
k^V-L ^ k^ ^ 1 (i.e. rj± ^ 1 and fc'^ 3> 1). In this limit the 
important effects are advection by the fluid displacement 
and magnetic diffusion. Generation of from Br by the 
Keplerian shear flow is negligible in this limit, but Coriolis 
and centripetal acceleration still play a crucial role through 
the dynamics of the fluid which enters via the appearance of 
the off-diagonal matrix elements in the 2nd term in eq ( |C1[ ). 

This short- wavelength, low-frequency limit corresponds 
to the cyclotron mode of the magnetised fluid, which has 
frequency 



eB^ fh 

r/nl miC p 



(C2) 



where the 2nd form applies for a simple ion-electron-neutral 
plasma ( [Wardle fc Ngl|1999| [Pandey fc Wardle||2008] ) This 
mode is able to couple effectively to the Keplerian rotation as 
long as the sense of circular polarisation of the mode matches 
the epicyclic motion, i.e. as long as B^rju < 0. The other 
short-wavelength mode, the high-frequency whistler (lj = 
fc'^v^/wijj^lis unable to couple effectively to the rotation 
| Wardle|[l999t . 

This limit applies for fc — >■ oo in regions II and III of 
Fig. |5] where arbitrarily large wavenumbers are unstable. 
The dispersion relation in this case is obtained by letting 
fc — >■ cxj in (|B6|), and reduces to 



(C3) 



(1 + v^){t]1 ryl) -f- |s77H 4- + 1 = 0, 



i.e. a.(y) = with a given by eq (B7l. The lack of any fc 



dependence in this regime occurs because both the magnetic 
diffusion and the magnetic stresses on the fluid (which are 
responsible for the fluid displacement) scale as k? . 



'usive limit Finally, the 2nd term in ( B2 1 may be 
neglected when fc'^r/x ~ 1 3> fc^ (i.e. r]± ^ 1, fc^ ~ l/ri±). In 
this case instability relies on the keplerian shear flow gen- 
erating B^ from Br, and the tendency of Hall diffusion to 
convert B^ back into Br- This brings the potential destabil- 
ising effect of Hall diffusion in shear flows to the fore and 
shows that it is quite independent of rotational effects - 
i.e. the Coriolis and centripetal acceleration - that drive the 
MRI |Kunz|l2008j). Here our simplified geometry and large 
diffusion limit allow us to find simple analytic expressions 
for the growth rate as a function of fc for arbitrary diffusivity. 
Eqs (Bll and (B2l show that \Sv\ <C IAjbI, so that the 



Sv term can be neglected in the linearised induction equation 



(B3l, yielding: 



u 

3 
2 



SB 

B 



VP svu 
-sr/H VP 



SB 

l3 



(C4) 



with Sv still given by eq (Bl I. In this limit the field evolves 



in response to shear and diffusion, without significant feed- 
back from the perturbations that it induces in the fiuid fiow. 
The 3/2-bearing term (= —d\nQ,/ d\nr) is the only mani- 
festation of the Keplerian rotation law and cylindrical ge- 
ometry in this equatiorQ To generalise this expression we 
replace |n by a characteristic shear frequency v' that can be 
thought of either as accommodating different rotation laws 
in cylindrical geometry, or represents plane-parallel shear in 
a cartesian geometry. In the latter case r and (j) are mapped 
to —X and — y with a local velocity field v — v' xy; the z- 
axis is unchanged. Without loss of generality we may assume 
that v' > 0, and as before, we write B — sBz with s = ±1. 



Upon replacing the 4 term in (C4l by v' , we obtain the 



dispersion relation in dimensional form: 



(C5) 



This immediately shows that instability requires both v and 
sr^H to be nonzero and have the same sign, otherwise the 2nd 
term exceeds the right hand side and there are no accept- 
able solutions. The dispersion relation is easily solved for 
the growth rate 



V = \J k"^ sr\Y{(v' — k^srjti) — fc' 



VP ■ 



(C6) 



and the equilibrium is unstable to wave numbers running 
from up to a cutoff 



fcc — 



svh V 



where, as usual, the perpendicular diffusivity is 



Vi 



vl + vl- 



(C7) 



(C8) 



To fi nd th e fastest growing mode we play the usual trick: 
write eq (C5l as a quadratic in fc^ and find vq, the value 



of V that gives zero discriminant. This means that there is 
only one corresponding wavenumber fcp and therefore the 
point (fco,i^o) lies on the peak of the curve u{k). The wave 
number is found from the quadratic, which is easily solved 



* In dimensionless form, ujjj = l/|r?Hl and the whistler mode has 
frequency ui = k"^ /\r]ii\. 



^ The velocity perturbations, however, are still afTectcd by ro- 
tation (see eq |B1[ | but do not themselves feed back on the field 
evolution. 
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when the discriminant vanishes. Fortuitously, this procedure 
yields simple expressions for the maximal growth rate, 

and corresponding wave number 



(Cll) 



where we have made use of the identity 
?7± - VP _ m 

TjH V± + VP 

to neatly avoid delicate subtractions when |77h| ^ Vp- Note 
that 

fo = kf,v± , (C12) 



which we can then use in eq ( B4 1 along with | »7h | , Vp ^ 1 



to show that the relationship between the perturbed field 
components in the fastest growing mode is 

SB^ = 5^5^ SBy . (C13) 

VP+V-L 

The MRI proceeds in the ideal MHD limit because shear 
creates from Br, BrB^ stresses on the fiuid cause it 
to move radially inwards and outwards, creating more Br 
available to be converted to B^ by the shear. In the diffusive 
limit, the radial component of B is not generated by the 
response of the fluid to magnetic stresses, but because Hall 
diffusion converts B^ to Br - For more general field and wave 
vector configurations, ambipolar diffusion may play a similar 
role in assisting or suppressing the MRI, albeit hindered by 
dissipation ( |Kunz fc BaibuslpOOi) |Desch|[2004l |Kunz|[2008l 
Pandey & Wardle, in prep). 

This paper has been typeset from a T^TjX/ M^rjX file prepared 
by the author. 
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